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OS ■ ABSTRACT 



in 



We use the rotation curves of a sample of dark matter dominated dwarf and low-surface brightness 
(LSB) late-type galaxies to study their radial mass distributions. We find that the shape of the rotation 
curves is remarkably similar for all (both dwarf and LSB) galaxies in the sample, suggesting a self-similar 
density distribution of their dark matter (DM) halos. This shape can be reproduced well by a density 
profile with a shallow central cusp [p{r) cx 1/r 7 , 7 a 0.2 — 0.4] corresponding to a steeply rising velocity 
curve [v(r) ocr 9 , g ps 0.9 — 0.8]. We further show that the observed shapes of the rotation curves are well 
matched by the average density profiles of dark matter halos formed in very high resolution simulations 
, of the standard cold dark matter model (CDM), the low-density CDM model with cosmological constant 

■ (ACDM), and the cold+hot dark matter model with two types of neutrino (CHDM). This is surprising 
[ in light of several previous studies, which suggested that the structure of simulated dark matter halos is 

C . inconsistent with the dynamics of dwarf galaxies. We discuss possible explanations for this discrepancy 

and show that it is most likely due to the systematic differences at small radii between the analytic model 
proposed by Navarro, Frenk, & White, with 7nfw = 1, an d the actual central density profiles of the dark 
matter halos. We also show that the mass distributions in the hierarchically formed halos are on average 
consistent with the shape of rotation curves of dark matter dominated galaxies. However, the scatter of 
the individual profiles around the average is substantial. Finally, we show that the dark matter halos 

■ in our hierarchical simulations and the real galaxies in our sample exhibit a similar decrease in their 
Oh| characteristic densities with increasing characteristic radial scales and show increase in their maximum 
q . rotation velocities with increase in the radii at which their maximum velocities occur. 

Subject headings: cosmology: theory - dark matter: halos — galaxies: kinematics and dynamics — 
c/3 , galaxies: structure 

1. INTRODUCTION tionless hierarchical collapse in CDM- like models was re- 

The amount of luminous matter (stars and gas) in many /S 1 ^ S !' ll( ^i e( ^„ US l I l 1 ^ resolution N body simulations 

. , , . , 1 • • V ffi • Z. J 1 • (Dubmski & Carlberg 1991: Navarro, Frenk, & White 

spiral and irregular galaxies is not sufficient to explain !„„. ir .„_ , n TT -, TTT i , , T ^ TTT „_ N „, , , 

=3 . X i-+ j j if t +v, • + +• (T>n \ 1996, 1997, hereafter NFW96 and NFW97). The halo 

the amplitude and shape ol their rotation curves RCs). , ., ' , r . . , , ' x , ., 

' r™. . „ . , , j c ,/ density profiles were found to be cuspy (corclcss) and well 

1 his discrepancy is usually interpreted as evidence for the „,, , , r „ , m /, T r n ,i m /.\ 

r . , , j , / mrt , , j fatted by the following two-parameter profile (Nb W96): 

presence of an extended dark matter (DM) halo surround- J & r r \ / 

ing the visible regions of galaxies (e.g., Casertano & van p s 

Gorkom 1991; Persic, Salucci, & Stel 1996, and references Pnfw(^) = — — — ; — 3. (1) 



o 

On 



therein). The extent of the dark matter halos, estimated 



(r/r s ) (1 + r/r s ) 



using satellite dynamics, is ~ 0.2 — 0.5h 1 Mpc (Zant- ml , . ,. , , ,. 

1 htl', TnnA m • 1 1 nn^ n -I , , 1 1 he characteristic density, p s . and radius, r s , are sensitive 

sky & White 1994: Cangnan ct al. 1997; Zantsky et al. , , . , . , 'Tf , '. , , , 

iri l„N TT i.i 1 • 1 j. -i 1- r n i 1 to the epoch of halo formation and are tightly correlated 

1997 . However, the dynamical contribution of the dark .,, ,, , , . . , /-kj^-kt nn rnf c ^ 

' ■ m • . with the halo vinal mass (NFW 96,97 . Therefore, the re- 
matter can be substantial even 111 the very inner regions ,, . , ,. v ' , J . 

r , . , , ,. , ... r f j. suits of tfaese simulations suggest a corclcss and sell-similar 

ol galaxies: the observed rotation velocities of some dwarf . tJ rr ,„, , .,. .. ... . . 

, , ,. , • 1 . , T nT1 , , . . , ,, density structure of DM halos, with the vinal mass being 

and low-surface brightness (LSB) galaxies imply that DM . , ,. 

, ■ . c )■ f . „ri»\ rj the single scaling parameter, 

constitutes a dominant fraction (up to ~ 95%) ol dynami- ° 5.i r 1 • 1 1 

, .,, . , \ . , '. ^.fi 1 1 he structure of the inner regions of galactic halos was 

cal mass within the last measured point of their RCs (e.g., ,. , , „, -n • 1 ?ioo^ j at /mn^\ 

r, . „ t-i 1AOO ,. K „ . o n studied by Flores & Primack (1994) and Moore (1994), 

Cangnan & Freeman 1988; Martimbeau, Cangnan, & Roy , , , . , , ,. .. v ' K . ' l 

-.„„. , T-x-i i „ „ , lrin -N m, j 7 „ , who used high-resolution rotation curve measurements of 

1994; de Blok & McGaugh 1997). Ihese dark matter aom- 111 1 • . 1 1 r 1 mi 

, , 1 ■ rr • c i- 7 ■ several dark matter dominated dwarf galaxies. Ihc ccn- 

mated galaxies otter a unique opportunity for probing ai- , , , . , . r , , 

,, ,f 1 . , n „, , , . 1 u Ii tral density distributions m these galaxies were lound to 

rectly the density structure of DM halos which can be then , . . •', , . , , , , 1 , , , . 

, .., ,. .. j. ., .- 1 ii be inconsistent with the singular or cx 1/r behavior 

compared with predictions of theoretical models. , , , . ™ ' ',. r , , 

rp, 1 . -i 1 . rm , i , r j - j . predicted by equation (1). the scaling properties ol the 

1 he detailed structure of DM halos formed via dissipa- , , , , v / , , „ , °, , c , 

observed halos were analyzed by Burkert (1995, hereafter 

1 We assume that the present-day value of the Hubble constant is Ho = lOO/i km/s/Mpc. 
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B95), who pointed out that shapes of the density profiles 
of four dwarf galaxies analyzed by Moore (1994) are re- 
markably similar and are well fitted by the following phe- 
nomenological density profile: 

PB ^ = (l + r/r b ) P [l + (r/r b f] ' (2) 

Parameters p b and r b were found to be strongly correlated, 
in qualitative agreement with the predictions of hierarchi- 
cal models (B95). 

In this paper we study the observed density structure in 
a sample of dark matter dominated galaxies inferred from 
their rotation curves. Particularly, we test two predictions 
of previous simulations of hierarchical halo formation: (1) 
cuspy central density distribution and (2) self-similarity of 
the halo density structure. We then use results of high- 
resolution A-body simulations to compare the observed ro- 
tation curves with circular velocity profiles of dark matter 
halos formed in different structure formation models. 

2. DWARF AND LSB GALAXIES 

2.1. The sample 

We have compiled a sample of 10 dwarf and 7 LSB 
galaxies with measured rotation curves and published mass 
models for stellar, gas, and halo components (see Table 1). 
The dwarf galaxies were selected from different sources, 
whereas all 7 LSB galaxies were selected from the sample 
of de Blok, McGaugh, & van der Hulst (1996) (see Table 

In our sample we included only those galaxies in which 
the dark matter component was shown to constitute ;> 
85% of the total mass inside the last measured point of 
the rotation curve (in most cases with the maximum disk 
assumption). It is important to note that distances to all 
of the dwarf, and some of the nearby LSB, galaxies are 
quite uncertain (in some cases by a factor of two). While 
rotational velocity is a directly observable quantity, the 
physical scale of rotation curves must be computed from 
the angular scale using distance. Thus, any uncertainty 
in the distance propagates into uncertainty in the physical 
scale. This fact should be kept in mind when one makes 
a one-to-one comparison of observed and modeled rota- 
tion curves. The distances to the galaxies adopted in our 
analysis are listed in Table 1. We have adopted the best 
estimate of distance from the original paper, when it was 
available, or the distance quoted in Tully (1988). We have 
also included in Table 1 the best fit values of mass model 
parameters described in §2.2. 

2.2. Rotation curve analysis 

Analysis of the dark matter distribution is difficult 
in most galaxies due to ambiguities in the estimates of 
the stellar mass-to-light (M/L) ratios and the resulting 
dynamical contribution of the stellar component to the 
observed rotation velocities (e.g., de Blok & McGaugh 
1997; Bottema 1997; Courteau & Rix 1997; and references 
therein). However, the rotation curves of the galaxies in 
our sample are mostly determined by dark matter on scales 
> 1 kpc: the contribution from the gas and stars is neg- 
ligible and residuals between the observed rotation curves 

2 Vnfw(?*) Whir/r for large r, which has an approximate slope 



and contribution of DM are at the level of the observa- 
tional scatter of rotational velocities. Therefore, the dark 
matter distribution models can be directly fitted to the ro- 
tation curves of these galaxies without uncertain assump- 
tions about M/L ratios. 

B95 showed that the density distribution, ps(r), de- 
scribed by equation (2) fits the data very well over the 
entire observed range of scales. At large radii this profile 
falls off as p(r) oc r~ 3 , in accord with simulations of the 
CDM model (e.g., NFW96). However, the change of loga- 
rithmic slope at r <~ tq predicted by pb(t), equation (2), is 
faster than the change predicted by pnfwW, equation (1). 
Moreover, ps{r) has a flat core at small radii (r <C n>), in 
disagreement with the r _1 central cusp of pnfwM- For 
dwarf galaxies, the scale at which the density distribution 
is expected to become flat is quite small (<J 1 kpc) and is 
in fact below the current observational resolution. From a 
theoretical point of view the existence of a core is difficult 
to understand because hierarchical formation of halos is 
much more likely to result in cuspy central density distri- 
butions (Syer & White 1997). Therefore, we will consider 
the broader family of density profiles (Zhao 1996): 



PW (r/r )^l + (r/ro)«p-^)/«- K > 

Note that p(r <C ro) oc r~ 7 , p(r ro) oc r~^, and a 
characterizes the sharpness of the change in logarithmic 
slope. This family includes both cuspy profiles of the type 
proposed by NFW96 (a, (3, 7) = (1,3, 1) and the so-called 
modified isothermal profile (a,/3, 7) = (2,2,0), which is 
the most widely used model for the halo density distri- 
bution in analyses of observed rotation curves. It is also 
convenient to make direct fits with an analytic model sim- 
ilar to (3) for the velocity profile: 

V(r)=V t (4) 

{) *[l + (r/r 4 ) a ](9+b)/a' W 

where rt and Vt are the effective "turnover" radius and ve- 
locity and a parameterizes the sharpness of the turnover. 
The limiting behaviors are V(r > r ( ) a 1/r 6 and V(r <C 
r t ) oc r 9 . The peak of the velocity profile (4) occurs at 
the radius r max = nig/b) 1 ^, and V max = V(r max ) = 
V t {g/b) g / a [l + .g/6]-( fc +f)/ a . 

The existing rotation curve measurements, due to their 
finite resolution and extent, cannot be used to constrain all 
five parameters of the profiles (3) and (4). The inner radii 
of the observational rotation curves are well probed, so 
we can make a meaningful comparison with fitting func- 
tions having different inner density profiles 7. However, 
most of the galaxies in our sample have rotation curves 
that have not begun to decline at the outermost mea- 
sured point and thus have very little information about 
the asymptotic slopes /?, b. Given the uncertainties, we fix 
the outer logarithmic slope to the value suggested by the 
models 2 (1) and (2): f3 = 3, b = 0.34. Note, however, that 
as was noted by Burkert (1995) this value is also favored 
(to (3 = 2 of the pseudo-isothermal profile) by observed 
highest quality rotation curves of dwarf galaxies. For the 
same reason, the sharpness of the turnover, a or a, are not 

b ~ 0.34 for values of r near a typical virial radius. 
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Table 1 

The sample of dwarf and LSB galaxies 
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^max 
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' max 
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Galaxy 


M B 


7 — 1 1 

a kpc 


1 — 1 

km s 


10 s /i J M Q kpc 13 


1 — 1 1 

a kpc 


1 — 1 

km s 


i — 1 i 

a kpc 


i — 1 

km s 


/i _1 Mpc 


Reference 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 












Dwarf 












DDO 154 


-13.8 


2.0 


38 


0.33 


2.9 


76 


5.6 


47 


3.0 


1,2 


DDO 170 


-15.2 


4.4 


52 


0.13 


6.5 


105 


12.4 


64 


14 


3 


NGC 2915 


-16.8 


1.1 


69 


3.77 


1.7 


140 


3.2 


86 


2.5 


4 


IC 2574 


-15.7 


5.1 


65 


0.15 


8.2 


138 


15.6 


85 


2.0 


5 


NGC 5585 


-17.5 


2.3 


73 


0.93 


3.4 


148 


6.5 


91 


4.7 


6 


DDO 236 


-16.8 


4.4 


59 


0.16 


6.8 


122 


13.1 


75 


1.4 


7 


DDO 7 


-17.7 


3.8 


87 


0.49 


5.5 


176 


10.5 


108 


18 


8 


DDO 10 


-16.3 


3.5 


54 


0.23 


5.3 


112 


10.0 


69 


7.8 


8 


UGC 2684 


-13.7 


1.6 


41 


0.59 


2.6 


86 


4.9 


53 


4.1 


8 


DDO 34 


-15.7 


1.5 


54 


1.29 


2.1 
LSB 


109 


4.0 


67 


5.9 


8 


F568-1 


-17.5 


3.8 


97 


0.61 


5.5 


197 


10.5 


121 


64 


9 


F568-3 


-17.7 


6.1 


96 


0.23 


9.1 


196 


17.4 


121 


58 


9 


F568-V1 


-17.3 


4.9 


96 


0.36 


7.2 


194 


13.7 


119 


60 


9 


F571-8 


-17.0 


6.8 


121 


0.30 


10.1 


248 


19.4 


153 


36 


9 


F574-1 


-17.8 


8.8 


95 


0.11 


13.6 


197 


25.9 


121 


72 


9 


F583-1 


-15.9 


4.6 


71 


0.22 


6.7 


143 


12.8 


88 


24 


9 


F583-4 


-16.3 


6.0 


64 


0.11 


9.3 


133 


17.8 


82 


37 


9 



NOTES.- Col. (2) Mb, blue absolute magnitude; col. (3) best fit ro(see eq. [3]; the fitting procedure is described in §2.2); 
col. (4) Vb = V(r ); col. (5) best fit p n (see eq.[3]); col. (6) best fit r t (see eq.[4]); col. (7) best fit V t (see eq.[4]); col. (8) r max ; 
col. (9) V m ax = V(r max ); col. (10) distance to galaxy adopted in this study; 

REFERENCES.- (1) Carignan & Freeman 1988; (2) Carignan & Beaulieu 1989; (3) Lake et al. 1990; (4) Meurer et al. 
1996; (5) Martimbeau et al. 1994; (6) Cote et al. 1991; (7) Jobin M. & Carignan C. 1990; (8) van Zee et al. 1997; (9) de 
Blok et al. 1996. 
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constrained for all of the galaxies. However, a fair number 
of galaxies in the sample do show the turnover and thus 
can be used to constrain a. The plausible value of the 
parameter a = 2 was determined using rotation curves of 
these galaxies. We generalize this value to all of the galax- 
ies (which in no way contradicts the data, but is not, of 
course, a strict procedure) and thus currently we can only 
talk about a plausible range of a values as a "universal 
fit" (if any such universal value exist at all). For example, 
our results will not change drastically if we use a = 1.5 
instead of a = 2. However, a = 1 gives a poorer fit to the 
data. 

We fix the parameter 7 to 0.2: the value which best fits 
most of the observed rotation curves. 

The corresponding best-fit slopes of the profile (4) are 
{a,b,g) = (1.50,0.34,0.9). Note that g = 1 -7/2. With 
parameters a, (3, and 7 (a, b, g) fixed, we fitted the data for 
the remaining free parameters of the profile (3): po and ro 
(Vt and r t in eq.[4]). Our fits thus have the same number of 
free parameters as do profiles (1) and (2). Note that while 
the particular set of the parameters (a, (3, 7)= (2,3,0.2) 
used in the paper didn't result from a strict fitting proce- 
dure, it was motivated by all possible constraints of avail- 
able data. The only theoretically suggested value is that of 
(3 but it also seems to be favored by data (Burkert 1995). 
Hopefully, as new RC observations come along, they can 
be used to pinpoint the parameters a and (3 with a better 
accuracy. 




Fig. 1. — Rotation curves of (a) dwarf and (b) LSB galaxies 
(symbols) normalized to the best fit values of ro and rota- 
tional velocities «o at ro predicted by density model (eq.[3]). 
The solid line on both panels shows the analytic rotation curve 
corresponding to the density profile (eq.[3]) with (a,/3, 7) = 
(2,3,0.2). The rotation curves for different dwarf and LSB 
galaxies have virtually identical shapes, which is very well 
matched over the entire observed range of scales by the ana- 
lytic model. Note, that the RC of NGC 2915 extends outside 
the scale of the plot: the outer part of this RC can be seen in 



Figure 1 shows rotation curves of dwarf (a) and LSB 
(b) galaxies normalized to their best fit values of ro and 
to the rotational velocities vq at ro, predicted by analytic 
profile (3). The best fit values of r , po, Vo = V(r ), r t , 
and Vt are given in Table 1 for each galaxy in the sam- 
ple. The formal errors of each of these values are less than 
about <~ 2 — 5%. Figure la shows that all of the dwarf 
galaxies have rotation curves of virtually identical shape 3 
with a remarkably small scatter. The rotation curves of 
the two dwarf galaxies, DD0154 and NGC2915, cannot be 
described by a smooth density distribution model in their 
outer parts. The RC of DD0154 shows a decrease in ro- 
tational velocity in the three outermost observed points. 
Conversely, the RC of NGC2915 has a sharp upturn at 
;> 5ft. -1 kpc (or r/ro ;> 4.5). This upturn can be seen in 
Figure 2. The explanation of this peculiar behavior is not 
clear (see, however, Burkert & Silk 1997), but it is obvious 
that it cannot be explained by any smooth model for the 
mass distribution. Note, however, that apart from the pe- 
culiar outer regions, the rotation curves of both DD0154 
and NGC2915 have the same shape as the rest of the galax- 
ies. 

The shape of the galaxies' rotation curves is well 
matched by the rotation curve corresponding to density 
profile (3) with (a, (3, 7) = (2, 3, 0.2) or correspondingly to 
RC (4) with (a, b, g) = (1.5, 0.34, 0.9). This result is in per- 
fect agreement with Burkert (1995) who showed similarity 
of rotation curves for four dwarf galaxies (two of which, 
DD0154 and DDO170, were included in our sample). As 
was mentioned above, the similar fit by pb(t) (eq. [2]) 
proposed by Burkert (1995) is equally good. Note, how- 
ever, that our profile docs not have any flat core, whereas 
Pb(t) predicts such a core at r C rj,. The fact that both 
profiles fit the data equally well is easy to understand if we 
notice that ps (r) predicts a flat density distribution at the 
scales well below the observational resolution (<; 1 kpc). 
Thus, pb (r) and profile (3) can be virtually identical in the 
range of scales resolved in observations and thus provide 
an equally good fit to the data. 

Figure lb shows that the rotation curves of dark mat- 
ter dominated LSB galaxies are also well described by 
the same analytic density profile. The larger amplitude 
of scatter in the case of LSB galaxies can be explained 
by the larger observational errors associated with a given 
point of a rotation curve and thus most likely reflects ob- 
servational uncertainties rather than intrinsic scatter of 
the halo properties. Most of the LSB galaxies in our sam- 
ple are located at considerably larger distances than dwarf 
galaxies. Therefore, the dwarf galaxies have been observed 
with considerably higher resolution and smaller observa- 
tional errors than LSB galaxies. The estimated errors are 
typically 10 - 20% (de Blok et al. 1996), especially in the 
inner regions of galaxies (<; 10 kpc). 

One important issue is whether subtraction of baryon 
component (stars and gas) in the galaxies from our sample 
will affect results of the rotation curve analysis presented 
above. As we mentioned in §2.1, the combined contribu- 
tion of stars and gas is <J15% for all of our galaxies (<^10% 
in most cases). It is clear that ideally one has to subtract 



Figure 2. 

By the shape of a rotation curve we mean its particular functional form. For example, the shape of the rotation curve described by equation 
(4) is x 9 /(l + x a )( g+b ' ) / a (where x = r/rt). By saying that the RC shape is similar for all our galaxies, we mean that all their rotation curves 
can be described by this functional form with fixed values of parameters a, b, and g. 
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contributions of both gas and stars from the observed RC 
in order to get the mass distribution of the dark matter to 
an accuracy of better than 10%. However, it is well known 
that this is not an easy thing to do. For stars, the exact 
mass-to-light ratio is not known and we cannot convert vis- 
ible flux into the stellar mass without making additional 
assumptions (hence, "maximum disk controversy" ) . In the 
case of gas, we know exactly how to convert the 21-cm flux 
into the mass of gas (no mass-to-light uncertainty). This 
conversion, however, relies on other assumptions which can 
easily lead to uncertainties as high as 10% in the contri- 
bution of gas. For example, we need to know distance to 
the galaxy in order to make this conversion. We also need 
to have a reliable way to estimate the profile of molecu- 
lar gas to recover the dark matter density profile. The 
distances to the dwarf galaxies in our sample are very un- 
certain (often by a factor of two or more) and so is the 
conversion. It is not clear whether it makes any sense to 
subtract the gas with uncertainties of its contribution this 
big. After all, the motivation for use of dark matter domi- 
nated galaxies for this kind of analysis is to avoid dubious 
or uncertain correction procedures which are very unlikely 
to result in a better determination of the shape of DM 
density distribution. The distance uncertainty is not so 
severe for LSB galaxies and so conversion could, in princi- 
ple, have been done in this case. We have not done this for 
one simple reason: subtraction of the gas component could 
change any given point of rotation curve by at most 10% 
(in practice less than that). However, the errors associated 
with each point of RC are of the order of 10 — 20% (which 
combine both observational errors and assymetries in the 
rotation curves between receding and approaching sides) 
and it seems unlikely that correction due to gas subtrac- 
tion would improve or systematically change the answer 
(unless the observational errors are significantly overesti- 
mated). We have tested the effect of gas subtraction on 
the RC shapes by performing RC shape analysis for two 
galaxies (dwarf NGC 5585 and LSB F583 - 1) with and 
without subtraction of gas. The two galaxies have been se- 
lected to have a clearly visible turnover of the RC and to 
have a fairly high fraction of gas inside the last measured 
point of the rotation curve. This fraction is 8% for NGC 
5585 (Cote et al. 1991) and 5% for F583- 1 (de Blok et al. 
1996). The results of fitting the (a, 0, 7) = (2.0,3.0,0.2) 
model to the RC with and without gas subtraction result 
in very similar results: the difference in the best fit param- 
eters is <J 10% and corrected and uncorrected RCs have 
virtually identical shape. Note also that the dwarf galaxies 
that we used have on average a higher (or at least as high) 
fraction of gas (typically <~ 6 — 10%) than LSB galaxies 
(typically 3 — 7%). Therefore, if gas would introduce sys- 
tematic differences in the shape of RC, we could expect 
that the scatter in Figure 1 to be larger for dwarf galax- 
ies (some galaxies have much more gas than the others). 
Yet, the shape of rotation curves for dwarf galaxies is very 
uniform. 

We have repeated the fitting procedure described above 
using the analytic profiles (1) and (2). As was mentioned 
above, pb(t) results in a fit that is equally good to the fit 
by profile (3) shown in Figure 1. However, the analytic 
profile proposed by NFW failed to produce a reasonable 
fit to the data, as was indeed pointed out in NFW96 (see 
their fig. 12). The major difficulty with this profile, as was 



noted before by Flores & Primack 1994 and B95, is that 
the inner slope of the density distribution (7 = 1) is con- 
siderably steeper than implied by the rotation curves. The 
finite spatial extent of the data and incorrect inner slope 
of the profile (1) lead to implausible solutions of the x 2 ~ 
minimization procedure (the values of r s increase without 
convergence) . 

The observed similarity of the shapes of the rotation 
curves for seventeen different galaxies, selected solely on 
the basis of their dark matter content, and the remarkably 
small amount of scatter, implies that their matter distri- 
butions are self-similar in terms of the density structure. 
Of course, this includes both stellar and gaseous matter 
as well as DM. Both stellar and gaseous masses are uncer- 
tain because of uncertainties in the distance, mass-to-light 
ratio, and atomic-to-molccular gas ratio of each galaxy. 
To the extent that we can neglect the stellar and gaseous 
components (a subject that we intend to address in a sub- 
sequent paper), the self-similar rotation curves of these 
DM-dominatcd galaxies imply that they all have the same 
density structure. The question we now ask is whether the 
disagreement between this density structure and Pnfw(^) 
indicates a failure of CDM-type models? 

3. COMPARISON WITH THEORETICAL MODELS 

3.1. Numerical simulations 

We have used the new Adaptive Refinement Tree (ART) 
A-body code (see Kravtsov, Klypin, & Khokhlov 1997 
for details) to simulate the evolution of collisionless dark 
matter in the three cosmological structure formation mod- 
els: (1) standard cold dark matter model (CDM: f2o = L 
h = 0.5, <j% = 0.7); (2) a low-density CDM model with cos- 
mological constant (ACDM: A = 0.7, h = 0.7, a 8 = 1.0); 
and (3) a cold+hot dark matter model with two types 
of neutrino (CHDM; fl Q = 1 and Q. v = 0.2; h = 0.5; 
<T8 = 0.7; cf. Primack et al. 1995). Here f^o, ^a, and f2„ 
are the present-epoch values of the density of matter, vac- 
uum energy (as measured by the cosmological constant), 
and massive neutrinos, respectively. The rms fluctuation 
in spheres of radius 8/i _1 Mpc, as, was chosen to conform 
with the local abundance of galaxy clusters, for ACDM 
and CHDM models it is also in agreement with measure- 
ments of the cosmic microwave background anisotropy by 
the COBE satellite. The simulations followed trajecto- 
ries of 128 3 cold dark matter particles in a box of size of 
Lbox — 7.5/i -1 Mpc. In the CHDM simulation, two addi- 
tional equal-mass "massive neutrino" species were evolved, 
which brings the number of particles in the simulation to 
3 x 128 3 . To test for the possible effects of the finite box 
size, we have run an additional simulation of the ACDM 
model with the box size twice as large: L^ ox = \hli~ 1 
Mpc= 21.43 Mpc. We will denote the two ACDM simu- 
lations as ACDM7.5 and ACDM15 according to their box 
sizes. 

We have used a 256 3 uniform grid covering the entire 
computational volume and finer refinement meshes con- 
structed recursively and adaptively inside the high-density 
regions. The comoving cell size corresponding to a refine- 
ment level L is Axi = Axq/2 l , where Axo = Lbox/ '256 
is the size of the uniform grid cell (L = corresponds to 
the uniform grid). The increase of spatial resolution corre- 
sponding to each successive refinement level was accompa- 



6 



KRAVTSOV ET AL. 



nicd by the decrease of the integration time step by a factor 
of 2. The simulations were started at redshift Zi = 40 in 
the CDM and ACDM7.5 simulations and at z% = 30 in 
the CHDM and ACDM 15 simulations. Particle trajecto- 
ries were integrated with the time step of Aao — 0.0015 
on the zeroth-level uniform grid in the case of the CDM 
and ACDM runs and with Aa = 0.006 in the CHDM run. 
The time step on a refinement level L is Aa^ = Aao/2 L . 
The time step for the highest refinement level corresponds 
to ;> 40, 000 time steps over the Hubble time. Six refine- 
ment levels were introduced in the highest density regions 
corresponding to a cell size of Ax$ — QAQh^ 1 kpc. The 
dynamic range of the simulations is thus 256 x 2 6 = 16, 384. 
Note, that the resolution is constant in comoving coordi- 
nates which means that actual physical resolution is higher 
at earlier epochs (the halos were resolved with six refine- 
ment levels as early as z ~ 1). The refinement criterion 
was based on the local ovcrdensity of dark matter particles. 
Regions with overdensity higher than 5 = nth(L) 2 3 ( L+1 ) 
were refined to the refinement level L. Here, nth(L) is 
the threshold number of particles per mesh cell of level 
L estimated using the cloud-in-cell method (Hockney & 
Eastwood 1981). We have used values nth = 5 at all levels 
in the CDM and ACDM runs; for the CHDM run we have 
used nth = 10 at the levels L = 0, 1 and n t h = 5 for all 
of the higher levels. These values of threshold were sug- 
gested by results of the tests presented in Kravtsov et al. 
(1997); they ensure that refinements arc introduced only 
in the regions of high-particle density, where the two-body 
relaxation effects are not important. 

For the dark matter halos used in our analysis the spa- 
tial resolution is equal to ss 0.5 — 2h~ 1 kpc (corresponding 
to the 6th to 4th refinement levels). For each of the ana- 
lyzed halos, we have taken into account only those regions 
of the density and circular velocity profiles that correspond 
to scales at least twice as large as the formal resolution. 
The mass resolutions (particle mass) of our simulations are 
listed in Table 2, and are in the range of ~ (1 — 10) x 10 7 ft. _1 
Mq. Therefore a typical halo of mass <~ lO 11 ^ 1 Mq in 
our simulations contains several thousands of particles. 

These simulations are comparable in spatial and mass 
resolution, as well as in the box size, to those of NFW96,97. 
There is, however, a significant difference: our simulations 
are direct simulations of all DM halos in a given compu- 
tational volume, whereas NFW96,97 simulate with high 
resolution a handful of individual halos. The fact that we 
analyze a statistically large sample consisting of dozens of 
galaxy-size halos in each simulation allows us to make con- 
clusions about average halo properties and estimate the 
amount of cosmic scatter. A summary of the numerical 
simulations is given in Table 2. The parameters listed in 
this table are defined in the text above. 

3.2. Tests of numerical effects 

There are several effects which can affect the halo den- 
sity profiles at scales larger than some particular scale re- 
lated to this effect. We have tested the reliability of the 
simulated density and velocity profiles by comparing re- 
sults of the simulations with different resolutions and time 
steps. Specifically, the tests were used to determine the 
range of numerical parameters for which the convergence 
of density profiles was found at scales larger than two for- 
mal resolution elements (formal resolution is equal to the 



size of the refinement mesh cell). 

Tests presented in Kravtsov et al. (1997) show that 
the density profiles are not affected by the force reso- 
lution down to a scale of about one resolution element 
(a similar conclusion was reached by NFW96). To test 
the effects of the time step we have used a set of 64 3 - 
particlc simulations of the CDM model with parameters 
identical to those described in the previous section. These 
test simulations were started from identical initial con- 
ditions, but evolved with different time steps: Aao = 
0.006,0.003,0.0015,0.00075. Comparison of the density 
profiles for the same halos in these simulations shows that 
for halos of all masses, the profiles converge for runs with 
Aa <, 0.0015 (the value used in our CDM and ACDM 
simulations) at all scales, down to the resolution limit. 
We further use two 128 3 -particle simulations of the ACDM 
model with the box size of 15/i _1 Mpc and with time steps 
of Aao = 0.006 and Aao = 0.0015. The comparison shows 
that the most massive halos (virial mass M v i r > lO 13 /^ 1 
Mq) have systematically shallower central (r <; 10 — 20h~ 1 
kpc) density profiles in the Aao = 0.006 run as compared 
to the halos from the Aao = 0.0015 run. However, the dif- 
ference decreases with decreasing halo mass and for M V i r 
<J 5 x lO 12 ™ -1 Mq the density profiles from the two runs 
are identical within statistical noise at scales larger than 
one resolution clement. This mass dependence is due to 
the different accuracy of numerical integration in objects 
of different masses. The accuracy depends on the average 
displacement of particles during a single time step: for 
the integration to be accurate, the displacement should 
be <J 10 — 20% of the resolution element. Particles in- 
side more massive halos have considerably higher veloci- 
ties (v ;> 300 — 400 km/s) and thus average displacements 
that are larger than the displacements of particles inside 
small halos (v <J 200 km/s). In this study we focus on 
the mass distribution of the small halos (M <J 1 x 10 12 
Mq), for which the tests indicate convergence of the den- 
sity profiles for time steps Aao < 0.006. The time step 
of all simulations presented in this paper, except for the 
CHDM simulation, is four times smaller than the above 
value (sec Table 2). As an additional test, we have com- 
pared average RC shapes for CDM halos in the 7.5/i _1 
Mpc box simulation shown in Figure 2a and for halos in 
the same mass range (<; 1 x 10 12 M ) from an identical 
simulation (identical initial conditions and simulation pa- 
rameters) with time step Aao = 0.006. We have found 
that average RC shapes and the scatter in these two sim- 
ulations are indistinguishable. 

The mass resolution in our simulations (particle mass) 
is - (0.6 - 5) x Itfh- 1 Mq for L box = 7.5/j- 1 Mpc 
runs, and 1.3 x 10 8 /i _1 M© for the test L box — 15/i _1 
Mpc ACDM run (see Table 2). Therefore, halos of mass 
M vir = 10 12 /i _1 Mq and M vir = 10 n /i _1 M Q (the range 
of masses used in our comparison with the data) are re- 
solved with <~ 100,000 and <~ 10,000 particles, respec- 
tively For reference, there are J> 100 — 200 particles inside 
the innermost point (2 formal resolutions) of the rotation 
curve used in the fitting procedure described below. Com- 
parison of the average velocity profiles in the ACDM7.5 
and ACDM15 simulations (the latter has eight times worse 
mass resolution than the former) shows that there are no 
systematic differences between profiles in these two runs 
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Table 2 

Parameters of the Numerical Simulations 



Model 


a 


n A 




h 


C8 


zt 


Aa 


L box 


Particle mass 
















xlO -3 


h- 1 Mpc 


xl0 7 /i _1 M 


CDM 


1.0 


0.0 


0.0 


0.5 


0.7 


40 


1.5 


7.5 


5.56 


ACDM 7 5 


0.3 


0.7 


0.0 


0.7 


1.0 


40 


1.5 


7.5 


1.67 


ACDM15 


0.3 


0.7 


0.0 


0.7 


1.0 


30 


1.5 


15.0 


13.3 


CHDM 


1.0 


0.0 


0.2 


0.5 


0.7 


30 


6.0 


7.5 


4.44f 0.56 b 



a cold particles; h hot particles. 



(see Fig. 2). 

The force resolution can introduce errors in rotational 
velocities. To estimate this effect, we assume that the fi- 
nite force resolution results in a flat core (p =const) at 
scales smaller than the resolution element h r in an oth- 
erwise ideal NFW halo (Eq. [1]). This results in the ve- 
locity profile v so f t (r)/v s = x //(x)/x/(l), where f(x) = 
x 2 h /3{l+x h ) 2 +F(x)-F(x h ), F(x) = ln(l+x)-x/(l+x), 
x = r/r s , Xh = h r /r s , v = v(r s ), and r s is the scale-radius 
of the NFW profile (Eq. [1]). This profile can be com- 
pared with the velocity profile corresponding to Eq. (1): 
v N Fw{r)/v s = y/F(x)/xF(l). The error is - 18% at 
r w h r , and <J 5% at r ;> 2h r (see Fig. 5). Thus, the ve- 
locity profiles of simulated halos should not be significantly 
affected at scales r J> 2h ri which is where we perform the 
fit to analytic models. 

To test whether the box size of our simulations (7.5/i _1 
Mpc) is large enough not to miss all important tidal effects, 
we have compared the density and velocity profiles of halos 
formed in ACDM7.5 and ACDM15 simulations. We have 
not found any systematic differences between halo density 
profiles in these simulations. The average profiles of ha- 
los are identical within the statistical noise (see Figure 2). 
We have also used another indirect way of testing for the 
proper simulation of the tidal fields. Tidal torques from 
the surrounding large-scale structure presumably play a 
major role in the acquisition of the angular momentum, 
J = |J|, by the galaxy-size halos (Peebles 1969; Doroshke- 
vich 1970). Therefore, we can test if the tidal effects were 
simulated properly by comparing the so-called spin pa- 
rameter for the halos in our runs with previous results 
based on the larger-box simulations. The spin param- 
eter, A, of a halo is defined as A = ^E] 1 ' 2 /(GM 5 / 2 ), 
where J is the angular momentum of the halo, E is 
its total energy, and M is the halo virial mass. We 
have found that the distributions of A is very nearly log- 
normal 4 , P{\) = (l/AV277cr)exp (-ln 2 (A/A*)/2er 2 ) , with 
A* w 0.047, 0.045, and 0.048 and a w 0.66, 0.55, and 0.62, 
for the CDM, ACDM, and CHDM models respectively. 
Our results are in good agreement with previous studies 
(e.g., Barnes & Efstathiou 1987; Warren ct al. 1992; Cole 
& Lacey 1996; Thomas et al. 1997). We therefore con- 
clude that our simulations properly include all essential 



tidal effects. 



3.3. Results 



We have used a halo-finding algorithm described in 
Kravtsov et al. (1997) to identify halos in the simulations 
at z = 0. The algorithm identifies halos as local max- 
ima of mass inside a given radius. The exact center of a 
halo is found iteratively. We have run tests to ensure that 
we determine the halo center and resulting central density 
profile correctly. Only halos with more than 1000 parti- 
cles within their virial radius 5 , r vir , were taken from the 
full list. Also, to avoid effects of ongoing mergers we ex- 
cluded those halos which have close (r < r 2 oo) companions 
of mass more than half of the halo mass. 

The circular velocity profiles, v{r) — [GM(r) /r] 1 / 2 , 
were constructed by estimating the mass inside concen- 
tric spherical A = lh" 1 kpc shells around the halo center. 
To avoid contamination by gravitationally unbound back- 
ground particles, we iteratively remove all particles which 
have velocities relative to the velocity of the halo as a whole 
greater than escape velocity. The escape velocity, v e (r) , at 
a given distance r from the halo center is computed analyt- 
ically assuming that density distribution follows NFW pro- 
file (eq. 1): v e {r) = y/-2<j>(r) w 2.15 v max ^ln(l + 2x)/x, 
where </>(r) is halo gravitational potential at the distance 
T ~i Vmax (r-max) is the maximum rotation velocity of halo 
and x = r/r, max . This is a good approximation for a 
large range of r: NFW96 show that profile (1) approxi- 
mates the density and velocity profiles of halos reasonably 
well at scales ~ (0.01 — l)r V i r (see also Fig. 5). The maxi- 
mum rotation velocity v m ax and corresponding scale r max 
are found from a halo velocity profile at each iteration. 
The analysis shows that unbound particles affect at most 
the outer regions of halos 6 , r <~ (0.5 — l)r viri while in- 
ner regions are virtually unaffected (Klypin, Gottlober, & 
Kravtsov 1997). 

To compare the shape of the observed and simulated 
dark matter velocity curves, we have fitted the latter with 
the analytic model described by equation (3) . The param- 
eters a, /?, and 7 were fixed at the values used to fit dwarf 
and LSB galaxies - 2, 3, and 0.2 correspondingly, and we 
have fitted for the remaining two free parameters - ro and 
Po. It should be noted that the observed rotation curves of 



4 Note that the most probable value of this distribution is Apc^ = A* exp[— a 2 ] ~ 0.032 for the fit values quoted. 

5 According to predictions of the spherical top-hat collapse model we define the virial radius as radius of a sphere encompassing a mean 
overdensity of 200 for CDM and CHDM models, and 340 for ACDM model (Lahav et al. 1991; Kitayama & Suto 1996; Gross ct al. 1997). 

6 The case, for example, for a small halo located at the outskirts of a larger system, or for two halos passing close to each other. 
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most dwarf and LSB galaxies are measured only to radii of 
<J 10 — 30/i -1 kpc and often are still rising at the last mea- 
sured point. Therefore, the mass distribution in the outer 
parts of the galactic halos (and often maximum rotation 
velocity) is poorly constrained. To avoid any bias in the 
fitting procedure we considered only the inner 30/i _1 kpc 
of the simulated halos. We then normalized each rotation 
profile to its best fit values of r and rotational velocity v 
at r and computed the average of these normalized pro- 
files over all halos considered in each cosmological model 
(~ 50 — 60). In Figure 2 we compare the average nor- 
malized dark matter velocity profiles for halos formed in 
CDM, ACDM, and CHDM models, shown by solid lines, 
with corresponding profiles of the dwarf galaxies from our 
sample, shown with different symbols (the symbols are as 
in Fig.l). The average profile from the larger-box ACDM45 
simulation is shown with a dashed line. This profile does 
not extend to values of r/r which are as low as for the 
ACDM7.5 profile (due to worse spatial resolution). How- 
ever, for values of r/ro, where the two profiles overlap, 
they are indistinguishable. The dotted lines show the 2a 
envelope representing the scatter of individual halo profiles 
around the average. It should be noted that the scatter in 
the inner regions of the halo velocity profiles is substantial. 
This scatter possibly reflects physical differences between 
individual halos: our tests show that it is unlikely that the 
scatter can be attributed to the statistical noise associated 
with the finite mass resolution. The mass resolution of our 
simulations is very high (see Table 2): the number of dark 
matter particles inside the smallest scale, r min , of rota- 
tion curve used in the fitting procedure is J> 200 for large 
(~ 10 12 M©) halos and £ 100 for smaller (~ 10 11 M©) 
halos. 

Figure 2 shows that on average the velocity profiles of 
halos formed in hierarchical structure formation models 
and observed dark matter halos are in good agreement. 
It also shows that both cold dark matter halos and halos 
of dark matter dominated galaxies exhibit a certain self- 
similarity of the mass distribution in their inner regions. It 
was noted previously (e.g., B95, NFW96) that hierarchi- 
cal formation of the halos should also result in well-defined 
scaling properties of the mass distribution. It is thus inter- 
esting to compare the scaling properties of galaxy halos in 
our sample with those of the DM halos formed in the three 
hierarchical models studied in this paper. Figure 3 shows 
the plot of the best-fit parameters ro and po of the model 
density distribution (3) for the dwarf (solid circles) and 
LSB (open circles) galaxies together with corresponding 
parameters of DM halos formed in CDM (a), ACDM (b), 
and CHDM (c) simulations. As before, the values of the 
remaining parameters of the profile (3) were fixed to (a, 
(3, 7)=(2,3,0.2). For both galaxies and simulated halos, 
the parameters r and p are clearly correlated: the ha- 
los that are compact are systematically denser. DM halos 
in all models are fairly consistent with the observational 
points, except possibly for the CDM model that appears to 
form halos somewhat denser than observed. Note that the 
absence of halos at r <, 2ft, -1 kpc is due to our finite nu- 
merical resolution rather than the generic failure of these 
models to produce very compact halos. The characteristic 
density of the DM halos correlates strongly with halo mass 
in a way that reflects the mass dependence of the epoch of 



halo formation (NFW96): low-mass small halos collapse 
at systematically higher redshifts (when the universe was 
denser) and are therefore denser than the larger higher- 
mass halos. Thus, the correlation observed in Figure 3 is 
likely to reflect the different formation epochs of individual 
halos. 




X 

0.1 I— — ■ I 

0.1 1 10 

r/r 

Fig. 2. — Average normalized dark matter velocity profiles 
for halos formed in (top panel) CDM, (middle panel) ACDM, 
and (bottom panel) CHDM models with corresponding profiles 
of the dwarf galaxies from our sample. The dotted lines show 
the 2a envelope representing scatter of individual halo profiles 
around the average. It should be noted that although the ve- 
locity profiles of the hierarchically formed dark matter halos 
are on average consistent with the shape of observed rotation 
curves, the scatter in the inner regions of the halo velocity 
profiles is substantial. This scatter possibly reflects real phys- 
ical differences between individual halos. The average profile 
from the larger-box ACDM15 simulation (with 2 times worse 
spatial and 8 times worse mass resolutions) is shown with a 
dashed line in the middle panel. This profile does not extend 
to values of r/ro which are as low as for the ACDM7.5 pro- 
file (due to worse spatial resolution). However, for values of 
r/ro, where the two profiles overlap, they are indistinguish- 
able. This suggests that the shape is not affected by the finite 
size of the simulation box and mass resolution. The peculiar 
upturn in the rotation curve of NGC2915 is discussed in §2.2. 

A similar correlation can be observed in the r max — v max 
plane, shown in Figure 4 (values of r max and v max for each 
galaxy are given in Table 1). The maximum point in a 
galaxy's DM velocity profile and the corresponding radius 
is a nice set of physical parameters for comparison with 
simulations. Ideally, such a comparison would not force 
any pre-supposed fit to either the data or the the simu- 
lated profiles. Unfortunately, most of the galaxy rotation 
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curves in our sample do not extend to r large enough to 
explicitly define the maximum velocity. Therefore, we find 
Vmax by fitting the velocity profile V(r) in eq. (4) using 
parameters which produce a velocity curve equivalent to 
the p(r) fits discussed above: (a,b,g) = (1.5, 0.34, 0.9). 
After each galaxy is fit with this profile, we use the fit to 
determine its maxima. We have, in a sense, tested this 
procedure with a mock run with the simulated halos: we 
performed the same fit to the inner profiles of the simu- 
lated halos (r< 30h _1 kpc) and then compared this to the 
maxima determined using a smooth (all parameters free) 
fit over all radii of the halos. The values of v max found by 
these two methods are virtually identical. Although there 
is some scatter between the corresponding values of r max , 
there is no systematic difference between the two. The rms 
difference in values of r max (determined by the two fitting 
procedures) is smaller than the scatter from halo to halo 
at any particular v max . 

We see from Figure 4 the trend that larger r max cor- 
respond to larger v max . Such a trend was also seen in 
NFW96 (their Fig. 10). Note that when the h factor is 
scaled out of the r max axis, the three models lie very closely 
along each other in this plot. But since the observation- 
ally determined radii of the galaxies depend on h, the plots 
show the differences accordingly. 

Note that we have fitted 7 = 0.2 from the observed 
galaxies rather than a value derived from the simulations 
because of the larger scatter in the latter. However, if we 
fix a — 2 and (3 = 3 we can attempt to find the best-fitting 
ranges of 7 for each of our simulated models. The proce- 
dure for finding best-fit gamma values for the halos is not 
entirely straightforward since, depending on the the mass 
of the halo (and intrinsic scatter), the characteristic scale 
ro, may be as small as the resolution limit, and the inner- 
slope will be poorly constrained. That is, if we were to 
fix a = 2 and /3 = 3 and fit the value of 7 freely for each 
halo, the smaller halos (M <; 10 11 M Q ) would not offer 
any strong constraint, and the fitting routine could yield 
spurious or unphysical values for 7. What we have done 
instead is compared the "stacking" of halos for a range of 
7 values (0-1.5). We fit each halo using logarithmic radial 
binning out to the virial radius, and only fit radii which 
were greater than twice the formal resolution. This bin 
spacing emphasizes the inner profile shape. After fitting 
each halo for a particular value of 7 (with (a, (3) — (2,3)) 
we normalized each by ro and po and effectively plotted 
them on top of each other. The scatter in this stack is a 
nice way of determining how well the imposed profile shape 
fits a set of halos. This procedure also allows us to treat 
large halos (which, within resolution, probe r/r <C 1) 
and small halos equivalcntly. Using this criteria, we find 
that the ranges 7acdm ~ — 0.4, 7cdm ~ — 0.5, and 
7chdm ~ — 0.7 provide the best fits. 

4. DISCUSSION 

The most important conclusion evident from the results 
presented in the previous section is that there is no sta- 
tistically significant discrepancy between the shapes of ro- 
tation curves of simulated halos and rotation curve shapes 
of galaxies in our sample. This conclusion is somewhat 
surprising in light of the previous results (Flores & Pri- 
mack 94; Moore 1994; NFW96) that indicated a signifi- 
cant discrepancy between numerical simulations and ro- 



tation curve measurements. We see at least two possible 
explanations for the controversy. First, both Moore (1994) 
and NFW neglected the fact that distances to these galax- 
ies (and thus the physical scales of the rotation curves) 
are very poorly determined (with a typical error of ~ 50% 
or more; L. van Zee, private communication). With a 
distance uncertainty this large, it is hardly legitimate to 
make a raw comparison of the profiles at a specific physical 
distance scale. The comparison between rotation curves 
of dwarf galaxies and different analytic models was made 
by B95. The comparison was made, however, in units 
of r/r;, (see eq.(2)), in which case the uncertainty in dis- 
tance in both r and r& cancels out in the ratio (the same 
is true for our analysis that was done in units of r/r ). 
Indeed, the discrepancy between the analytic model (1) 
and the data (Fig.l in B95) was not as large as was found 
by Moore (1994). Note that the study by Flores & Pri- 
mack (1994) was also done using dimensionless scale units: 
r/biii, where r is the physical scale of the rotation curves 
and bni is the HI disk scale length. The key point is that 
the rotation curve shape is independent of the uncertain 
distance to the galaxy. 

The second possible source of discrepancy concerns the 
procedure followed to compare the numerical results with 
observations. When comparing to the RC data, we have 
used (a, f3, 7) = (2,3,0.2) fit (motivated by the observed 
galaxies' RCs; see Fig.l) only to determine r and vq for 
the inner halo RCs (i.e., r < SOh^ 1 kpc) so that these 
can be properly rescaled and compared to the dwarf RCs 
in Fig. 2. Other authors compare the data to a simple 
analytic fit to the entire halo profile (i.e., r <^ r V i r ), e.g. 
PNFw(f)- Therefore, when compared to the data, any 
deviations of actual halo profiles from analytic fit were ne- 
glected. Although the universality of the mass distribution 
in the DM halos is most likely real and reflects the self- 
similar nature of their formation, the associated scatter of 
the real profiles should not be neglected. Also, possible 
systematic deviations (especially in the inner, r <^ ro, re- 
gions) of the actual profiles from the analytic model (1) 
should be kept in mind. Our analysis shows that such de- 
viations do, in fact, exist. Figure 5a shows the velocity 
profiles of a sample of the DM halos in our CDM simula- 
tion normalized to their best fit values of the characteristic 
radius r s and rotational velocity at r s . Figure 5b shows 
residuals between the halo velocity profiles and the ana- 
lytic fit by the NFW profile. All profiles are shown down 
to their spatial resolution. Figure 5b shows that the halo 
rotational velocities at scales r/r s <; 0.5 are systematically 
lower than the rotational velocities predicted by the best 
fit NFW profile. Fig. 5 shows that pnfw{ t ) is a rather 
good fit to our halos for r ;> 0.03r„j r , but not for the 
smaller scales that are relevant to observed inner rotation 
curves. Figure 6 in NFW96 and Figure 4 in NFW97 show 
that similar deviations seem to exist in their simulations 
as well. It was suggested (J. Navarro, private commu- 
nication) that the inner density distribution may depend 
on the dynamical state of the halo: the most relaxed ha- 
los may have systematically steeper inner density profiles. 
We do not find such a trend for the halos analyzed in this 
paper. There does not seem to exist any correlation of the 
inner slope of the density profiles or the concentration pa- 
rameter, c = r vir /r s , with the dynamical state of the halo 
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Fig. 3. — Correlation of the best-fit parameters ro and po 
for the dwarf and LSB galaxies (solid and open circles, respec- 
tively) and for the dark matter halos (crosses) formed in (a) 
CDM, (b) ACDM, and (c) CHDM models. The correlation is 
consistent with the correlations of simulated DM halos: smaller 
halos are denser. 
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Fig. 4. — Correlation of the maximum velocity of the rotation 
curves, v max , and scale at which this maximum occurs, r m ax, 
for observed galaxies and simulated DM halos (the symbols are 
the same as in Fig. 3). 



quantified by the fractional difference between the center 
of mass inside the halo virial radius and the halo center 
(the density peak): d C M = \r pea k - r C M\/r V i r - However, 
if a weak correlation does exist, it could be lost due to the 
rather large errors in determining c (~ 30%) for small ha- 
los (M vir ~ 10 11 - lO 12 /^ 1 M ). Our analysis shows that 
although the whole ensemble of the halo profiles can be de- 
scribed reasonably well by a fixed set of the parameters a, 
0, and 7, the scatter in these parameters among the indi- 
vidual halos is substantial. Thus, for example, the average 
velocity profiles of halos shown in Figure 2 are very close to 
the observational points, whereas the upper 2cr-envelopes 
(the dotted lines) lie considerably higher in the inner re- 
gions of the rotation curves. The discrepancy between the 
analytic model of NFW and the data found by B95 is actu- 
ally within <J la from our average halo profiles. The scat- 
ter is not caused by poisson noise: the halo rotation curves 
contain ~> 100 particles inside the innermost bin used in 
the fitting. There are different possible sources of this 
scatter. The set of parameters (a, 0, 7)= (2,3,0.2), used 
to make a sensible comparison between data and models 
in the r/r — v/vo plane may not be the best description 
of the mass distribution in the simulated halos. However, 



our analysis shows that the true source is probably real 
differences between mass distributions of different halos. 
For example, when we allow parameter a vary freely with 
other parameters fixed to 7 = 0.2 and = 3, this results 
in a wide range of best fit values for a (<~ 0.7 — 2). Also, 
the density distribution in the inner regions varies sub- 
stantially from halo to halo: the density profiles of some 
halos are considerably steeper than profiles of the others. 
We can thus talk only about "approximate universality" of 
the profiles and some amount of scatter will be introduced 
with any fixed set of parameters. Regardless of the nature 
of this scatter, it should not be neglected when making 
comparisons with the data. And it is this scatter which 
makes us conclude thus it is premature to claim a dis- 
crepancy between the mass distribution of hierarchically 
formed halos and observed rotation curves of dark matter 
dominated galaxies. 
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Fig. 5. — (a) Velocity profiles of the DM halos in our CDM 
simulation normalized to their best fit values of the NFW 
profile: the characteristic radius r s (see eq.[l]) and the ro- 
tational velocity at r s (the actual values of these parameters 
for the halos used are in the range r s ~ 4 — 10/i - kpc and 
v s ~ 120 — 200/i -1 km/s depending on the halo's virial mass). 
The dashed line shows the RC predicted by the NFW an- 
alytic density profile (eq. [1]). (b) Residuals between the 
halo velocity profiles and the analytic fit by the NFW pro- 
file. All profiles are shown down to their spatial resolution. 
Note that although the NFW profile provides a reasonably 
good fit at 0.5 < r/r s < 20 (corresponding approximately to 
0.02 <J r/r vir < 1), the rotational velocities of halos at scales 
r/r 3 < 0.5 are systematically lower than the rotational veloc- 
ities predicted by the NFW profile. The long-dashed curve 
shows errors in rotation velocity due to the force softening, 
assuming the NFW profile (see §3.2). The errors due to soft- 
ening are < 5% at the scales that were used in fitting (x > 0.2, 
corresponding approximately to two resolution elements) and 
cannot account for ~ 5 — 10 times larger deviations from the 
NFW profile observed at these scales. 

Although it appears that we cannot reject any of the 
analyzed cosmological models, we note that the cores of 
DM dominated galaxies are potentially very useful probes 
of both the history of galaxy formation and the underly- 
ing cosmological model. The inner regions of DM halos 
are expected to be very sensitive to the merging history. 
The accretion of small dense satellites builds up a cuspy 
inner density profile after only few mergers (Syer & White 
1997), even if the initial density profile had a flat core. 
The slope 7 is expected to be a function of the merger 
rate and the slope of the perturbation spectrum on scales 
corresponding to the mass of a halo. This result can be 



used possibly to constrain models. According to the anal- 
ysis of Syer & White (1997), the steeper the slope of the 
power spectrum n (P(k) oc k n ), the shallower the central 
density profile: 7 w 3(3 + n)/(5 + n). If this result is 
qualitatively correct, the shallower slope of P(k) at clus- 
ter scales (n ~ — 2), as compared to the slope at galactic 
scales (n <~ —2.5), should result in steeper central den- 
sity profiles for cluster-size halos. The rotation curves of 
galaxies in our sample suggest a slope of 7 ~ 0.2 — 0.4. 
These values, according to the analysis of Syer & White, 
correspond to spectrum slopes of n ~ —2.9 and n w —2.7, 
Note that ticdm ~ -2.4, racdm ~ «chdm ~ -2.6 on 
the scales of the halos that we analyze. The slope 7 may 
potentially be a useful probe of the spectrum. Note that 
our simulations show a considerable spread in the slopes of 
the central density profiles. The physical processes which 
lead to differences in the central density profiles seen in 
our halos (e.g., environment, dynamical state, etc.) are 
not clear and will be investigated in a future study. 

Navarro, Eke, & Frenk (1996) suggested that the cores 
of dwarf galaxies may also be sensitive to possible past vi- 
olent star formation bursts in these galaxies. They showed 
that a sudden loss of most of a galaxy's gas, blown away in 
the course of a star formation burst, may lead to formation 
of a flat core, even if the initial density distribution was 
coreless. However, the parameters used to model such an 
event were somewhat extreme, and it is not clear whether 
this mechanism would work for most of the dwarf galax- 
ies. Moreover, a significant fine tuning would be required 
to explain the observed degree of self-similarity of the mass 
distribution in the dwarf galaxies. In this paper we showed 
that mass distribution in the LSB galaxies is, in fact, very 
similar to the mass distribution of the dwarf galaxies. This 
makes it even harder to explain such regularity with this 
mechanism alone, because many of the LSB galaxies are 
fairly massive systems which are unlikely to have lost much 
of their gas in starbursts. 

Recently, Burkert & Silk (1997) used an improved ro- 
tation curve measurement of the dwarf DD0154 to argue 
that discrepancy between the NFW profile and the ob- 
served mass distribution in this galaxy can be explained 
by a dark spheroid of baryons with mass several times the 
mass of the observed disk and comparable to the mass of 
the cold dark matter halo. This hypothesis must be tested 
using other galaxies and observations of the MACHOs in 
our galaxy that constitute the main observational evidence 
for the existence of such massive baryonic halos. The ob- 
served self-similarity of the mass distribution in the dark 
matter dominated galaxies which we analyze requires that 
the distribution of the baryonic dark matter be self-similar 
too and suggests a certain degree of "conspiracy" between 
baryonic DM and cold DM halos. The results presented 
in this paper show that the mass distribution in our hi- 
erarchically formed cold dark matter halos is consistent 
with the dynamics of the dark matter dominated galaxies 
and we think that it is somewhat premature to invoke an 
additional dark matter component. However, if further re- 
sults on the distribution of MACHOs shows that the mass 
contribution of these objects is substantial, the possibility 
of such a component must be taken into account. Note 
also that, as we discussed in §2.2, the rotation curves of 
two of the dwarf galaxies, DD0154 and NGC2915, show a 
peculiar behavior at large radii that cannot be described 
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by any smooth model of density distribution. 

5. CONCLUSIONS 
To summarize, our results are the following. 

• Rotation curves of DM dominated dwarf and LSB 
galaxies have a similar shape. This shape is incon- 
sistent with pNFw( r ) [1])) but is well fit by 
the coreless profile described by eq. (3) which has 
a shallower slope at small scales: p(r) oc 1/r 7 , 7 w 
0.2 — 0.4, corresponding to a steeply rising velocity 
curve [v(r) car 9 , g m 0.9 — 0.8]. 

• We find that on average the velocity profiles of the 
halos formed in the hierarchical structure forma- 
tion models analyzed in this paper (CDM, ACDM, 
and CHDM) and observed dark matter halos are in 
reasonably good agreement. We find a substantial 
amount of scatter in the central density profiles of in- 
dividual halos around the average. The physical pro- 
cesses which lead to differences in the central density 
profiles seen in our halos (e.g., environment, dynam- 
ical state, etc.) are not clear and will be the subject 
of a future study. 



• The inner (?* < 30ft- -1 kpc) average density profiles 
of DM halos in all of our simulations are well fit by 
model (3) with (a, f3, 7)=(2,3,0.2), and equivalently 
the rotation curves are well described by RC (4) with 
(a,b,g) — (1.5,0.34,0.9). The profiles systematically 
deviate from the NFW profile (1) at small scales. 

• We find that dark matter dominated dwarf and LSB 
galaxies show correlations between their character- 
istic density and radius consistent with the correla- 
tions of hierarchically formed DM halos: physically 
smaller halos are denser. We find a similar corre- 
lation between the maximum of the rotation curve, 
Vmax, and the corresponding radius r max . 
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